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AMSTMACT 


of •oditun and OH c oB ci ^ wi^otto la Ua 
f^/ 02 /H 2 naoMB by respoctivo si^aratod and &ow->powor lator 
flitoroBC«EMe teclmiqtiea bave permittod a aetailed eitaxxdaatioa the 
^<»oaiced flame chemistry of sodium In such oxygMi rich mec^. 
Previcms iater|Hreftaticms have been siu>wn to be lai^eiy incomplete or 
in error. 

The flame downstream inri^aies Indicate tl^t the amoimt of free 
sodium tracks the decay <d H>atom and as the flame radicals decay 
sodium becomes increasingly bound in a nuilecular form. A detailed 
kinetic model indicates that the sodium is distributed between NaOH and 
Na02 species. Concentratioas of NaO are very small andNaH negl^(ible 
The actual distribution is controlled by the state of equilibration o' the 
flames' basic free radicals. Na, Na02 and NaOH are all coupled to one 
another by fast reacticms which can rapidly inter convert one to another 
as flame conditions vary. Above aboiit 2000K, NaOH becomes dominant 
whereas Na 02 plays an increasingly important contribution at lower 
temperatures . 

The dissociation energy 0^(Na>0.) is established to be 35 + 5 kcal 

- 1 VC I y 

mol and the rate constant for the Na^iO^+M reaction of lx 10* cm” 
- 2-1 * 

molecule s , not incompatible with a recent flash photolysis measure- 
ment at lower temperatures . The rate constant for the reaction 

NaOH-<-H ^Na+HyO approximates to its maximum allowed value based 

^ -12 1/2 3 

on gas kinetic hard-sphere collisions, k - 1. 1 x 10 T cm molecule 
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INTROetFCnCN 


Al:Uu>agh it U dS£i^alt to a Hato* chaa^bt who haa aot atodtod 
aodiimi. thb llaaaa chaitoatry td thia «d»iqtdtOtta atom haa not ba^ aaUa- 
factorily charactariaad tmdmr many ccnnnMn flame ccHi^Uti<ma . Thia ia 
a conaaqnttkca ^ onr continuing toabil^ to mcmitor directly molacolar 
aodiom apaciaa in flamea . The anbjact of thia preamtaticm is the 
identificati.oa of the dtmiinant chemistry accounting for the observed 
aodinm disappearance prttfihss in a series of atmos|dieric pres sore, 
lean (oxygen rich) 8*^* flows of widely varying compo- 

sition and temperature. 

bx rich flames sodium remains essentially in elemental 

form.^' ^ Small amounts of NaOH formation can be expected via 

*^12 

Na + H^O T=i NaOH + H (12)* 

^-12 

which only becomes significant at temperatures above 2300K. hi lean 
112 / 02/^2 ihe decreased H-atom concentrations would favor 

si4(htly increased NaOH producticm. 

3 

Kaskan monitored Na and OH profiles in absorpticm in a series of 
lean, atmospheric pressure flames at 1400 to 1700K. He 

found reacticm (12) not balanced and unable to account for the measured 
decays. Kather, he found the initial N a decay rates proporticmal to the 
O 2 concentration and sv^gested the three-body reactimi 

‘'i 

Na +0^ +M Na02+M (1) 

k-1 


* Reaction numbering corresponds to ordering in the master list of 
reactions of Table n. 
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as Om danlmat Ms ozlistlon |»N»esss ta tiiisss Gmmmm . ISs sbmsw^ 

.«* ^ _2 

dbcsy rates eorrsspmd to a vatea of ab^^ 10 cm aureate 
s'^^ for k|. McEwan and te a s^s^o^t stody with, teiapsra* 

taros rate(iim from 1380 to 2200K also moaiton^ sediom ia slNiorp^oa 
bat otilissd tte Li/ldOM aa«Uu>d to fiMasare H>atOTi. Tbsy fooi^ also 
that reaettem (12) coold not dsscribs th« obsorvmiMa decay pre^es. 
However at temperatorss above 2000K they cot^luded that reactioo (1) 
was balanced aadMa02 was dtmkinant sodium oaddation product. At 
lower temperatures it a{^»eared to be too slow to become balanced. Fr<an 
the temperature dependence of a |dot of ^aO 2 )/ 0 ^a)(O 2 ) against they 
calculated a value of 65 kcal mol‘^ for DQ(Ma-02). This value was found 
to have a computational error and teter revised downward to 55.9 kcal 

mol* . ^ Frmn the initial Na decay rates in their flames McEwan and 

4 *33 6 

Phillips^ also determined values for k. of about 2 - 3 X 10 cm 

“2 “1 * 

molecule e in reas<mable agreement with Kaskan'^ values. 

ha an extension of Kaskan's earlier effort. Carabetta and Kaskan^ 
examined the sodium oxidation in a series of lean flames as a functi<m of 
pressure to evaluate the third body dependence of the Na disappearance 
rates. Over the temperature range from 1420 to 1500K at pressures from 
100 to 1520 torr they found the initial rates of decay of sodium to be pro* 
portional to both O 2 and M concentraticxis supporting reaction (1) as the 
dominant sodium removal process. 

7 

Recently. Husain and Plane measured Na oxidation rates in a heated 
cell at 724 and 844K. They generated atomic Na by flash photolysis of 
Nal and followed the Na decay in absorjXion in the presence of He, N 2 . 

and CO, baths. They found theNa removal rate described by reaction (1} 

“ -30 6 “2 -1 

with a rate constant of about 10 cm molecule s~ with N 2 the 

third body. This is 3 orders of magnitude faster than the reported 

measurements for k^ in the previous flame studies . The unambiguity of the 

static cell measurements questiems the assignment of the measured sodium 

decay rates in flames to reaction (1). Obviously some global overall 

process has been observed in flames that appears to behave as reaction (1). 
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orwonfliwu" 

Vary rtcwntly, Jwu«b aad J<pm‘ cianliMd tha eatalyaU of Hum 
nt/ficml r«cotn^i«tioa r*t«« hf ei sodiom in a of 

fttol'Heh ^^^i^**^** Tf>/f tsmA could ox^in Hh» sodtuai 

cotalysod H-atcnn <U«o|^^rai^o rotoc with rMCtloo (U) atom wi^ 


Ka-tOH-tht »• MaOHm 


-27 -1 

A ccmipater fit to their measoremrats yielded a value of 5 N 10 T 

cm^ molecule”^ e**^ for k.,. This correspcmds to 2.5 x 10*^^ cm^ 

-2-1 *■* 9 

molecule a at 2000K, unexpectedly large. Jensen then found that 

3 

using reacUcms (12) and (1?) could account approximately for i^kan's 
lean flame Na oxidaUcn rates without recourse to the inclusion NUO 2 > 
This, of course, eliminates the discrepancy existing between values for 
k^ measured in flames and in the static cell. Jensen concluded that for 
reaction (1) to have a negligible contribution to the flame observations, 
but still be consistoit with Husain and Plane's^ measurements, would 

require a value of about 40.6 kcal mol~^ for D^(Ha-0,), about 13 kcal 

-I ^ 5 

mol below the revised value cd Daugherty, Me Ewan and Phillips. 

However, this test by Jensen utilized the very limited data base provided 

by Kaskan which involves only two very similar flames . Our recent 

measurements on ten flames provides now the opportunity for a more 

extensive modeling of this system and a more critical test of Jensen's 

hypothesis . 
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Most of ths expsrimentsl system has be«ii described previously 

and will be very briefly discussed below. Flames were stabilised on a 

12 

2 cm dia . flat -flame ^dley-Sugdea burner ccmstructed of bundled 

hypodermic tubing . The burner is manifolded such that separate metered 

H 2 /O 2 /K 2 flows can be fed to an inner bundle. 1 cm dia. . and to i-.e 

outer bundle. An aerosol of aqueous sodium chloride was added to the 

3 -1 

central burner core by passage of a few cm sec of throt^h an ultra- 
13 

sonic nebuliser. Atomic sodium is rapidly generated in the flame 

reacticm seme and is introduced at a rate to produce concentrations of 
10 -3 

about 10 cm in the hot flame gases. 

The oj^cal arrangement is illustrated in Fig. 1. ANd-YAG (Quanta 
Ray DCR-PDLi-1) pumped tunable dye laser was used to monitor Na and 
OH in fluorescence. A calibrated brightness temperature tungsten ribbem 
lamp was used for sodium line reversal measurements of flame temperatures . 

The laser beam, about 1 mm dia^ , was passed through the flame and 
fluorescence sampled at right angles with an £/5.3 spherical mirror to 
match the monochromator aperture. The sampled beam was passed through 
a two mirror image rotator and into the entrance slit of an 0 . 35 m mono- 
chromator (McPherson 218). The detection optics image a section of the 
horizontal laser beam into the vertical monochromator slit with unity 
magnification. Thus, using a SO/xm slit that is 0.5 cm high, we detect an 
0.5 cm length of beam fluorescence, near flame center, with a flame height 
resolution of the order of 100 um. The burner wits mounted on a platform 
that was raised or lowered relative to the system optical axes to obtain 
time dependent profiles of properties above the flame front which is 
positioned close to the burner. 


Sodium fluorescence measurements were made by saturating the 
2 2 

sodium 589.0 nm ( P^^ 2 ) monitoring the closely 

coupled fluorescence. Under saturated conditions the 

fluorescence is insensitive to quenching and variations in Laser power. 

With the short (6 ns) laser pulse there is insufficient time for either laser 
11 14 

induced chemistry or ionization of the sodium. In this manner the 
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Figure 1. Experimental arrangement for laser induced f1u;>r::;;««utfe 
and line reversal temperature measurements. 
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relative fliM>reecence intmeities were proportional to the attmiic aodkum 
coneeatratioaa. Relative concentratton nMaearexn«[Kta were adeqo^e 
for thia atody. 

Abaolute OH cc»c«citrati«&a were meaaored in flnoreaceace oaing 
onr previoualy eatabliahed experimental condititmaf^' The 

OH(A^ ** ) (1* 0), R^(6) tranaition at 2S1. 14 nm waa excited and 

flaoreacence detected at 314.69 nm for the (1, 1) Q|(7) tranaition. The 
doubled output from the dye laser (using an R590 dye) was attenuated 
to produce fluorescence in the low power regime where fluorescence 
intensity is proportional to laser power. These relative intensity 
measurements were rendered quantitative by calibrating against OH 
fluorescence measurema3.ts in a high temperature fully equilibrated 
flame 0^702^2 ~ for which the OH concentration is readily 

calculated. 

The concentrations of the other flan^ radicals and trace flame con- 
stituents, ({]), (O), and (H 2 ) were calculated assuming balanced 
chemistry in reactions (18), ( 19 )> and (20) 


H+02==*» 

OH +0 

(18) 

0+«2 — 

OH +H 

(19) 

OH+H 2 = 

K 2 O + H 

(20) 


Kinetic .calculations confirmed the ability of these processes to remain 
balanced throughout the flames employed in the study. 

A flame znatrix was employed to cover a wide range of compositions 
and temperatures. Flames studied are listed in Table I. These encompass 
a 750K temperature range from about 1650 to 2400K. In the hot flame 
gases the two major species H 2 O and O 2 have ccmcentrations that vary by 
factors of x 1.2 and x24 fold, respectively, over this range of flames, 
hi all cases significant departures from equilibrium are apparent for the 
flame radicals with extensive decay occurring over the first few milli- 
seconds of the burnt gas region. 
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. Experimental flame matrix and some major flame properties. 
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MEASURED COfCQf TRATION PROFILES 

The ccmceatrationi of free eodium have been obtained as a function 
of distance (time) throughout the burnt gas regio&s of the ten oxygen rich 
flames . The fluorescence intensities have to be corrected solely for the 
differing flow volumes of each flame because the same rate of addition of 
sodium is n.\de to each. Resulting concentrations over the first 4 ms regi<m 
are indicated in Figure 2 for selected flames. These indicate the prtmounced 
decnys as observed previously. * The decay shapes vary from flame to 
flame. This behavior is xnari^dly different from that in fuel rich hydrogen 
or acetylene flames where the sodium remains almost totally in its free 
atomic state with only slight losses due to some NaOH formation and to 
diffusion. At these temperatures, for our purposes, thermal ionizaticm is 
negligible on this time scale. Molecular formation is quite evident, be> 
coming especially pronounced for certain flame compositions. Also, there 
is the suggestion that each flame starts in the reaction zone region from the 
same level of totally free atomic sodium which is then subject to its specific 
rate of depletion. 

Corresponding OH concentrations in these flames are illustrated in 
Figure 3. Based on our earlier developed technique for quantitative low- 
power laser induced fluorescence measurements, which have been 

shown to be valid for fuel rich flames, this present series of measurements 
appear similarly valid. 

Numerical kinetic calculations concerning the detailed balance among 
the minor flame species in oxygen rich flames, namely for OH, H, O and H^, 
indicate, as under fuel-rich conditions, that there is sufficient time for these 
radicals to become interrelated in the burnt gas regions. By considering 
reactions (18), (19), and (20) it can be shown that, in the present time frame, 
changes in the concentration of one radical will immediately reflect through 
to the others. Consequently, from measured OH concentrations it is possible 
to calculate the corresponding amounts of H, O, and from the appropriate 
equilibrium constants and concentrations of H2O and O^. These latter two 
major species are largely unperturbed from their equilibrium concentration 
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[Na], (relative scale) 



DOWNSTREAM TIME, (ms) 


Figure 2. Meaiured atomic sodium concentrations in the burnt gases 
of several oxygen-rich hvdrogen flames. Total sodium 
input approxiniates to 1010 atoms cm'^. 
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Figure 3. Sxperimentmlly measured OH ccmcentrations in the burnt 
gases of several oxygen-rich hydrogen flames. 
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▼mlMi and «o can ba ealcolatad by an aqalUbrlwn coda at tha xnaaanrad 
flama tamparatora . fa tUs cray» tba lolloadng aaq^raaaioaa can ba darfrad: 

1«J • K,gK„i6jj [k,6] • 1°1 * Tq;-n^-a] • 

r , 

[«2l • k,gK- I6gj • 

Tha relationship between OH and H is markedly different from fuel-rich 
behavior where [OK] a [H] a [Oj^^^o; [ 02 ]^^^# This is apparent in the shapes 
and magnitadeu of tha [H] decays in Fignre 4 calculated from the corraspimd- 
ing measured OH concentrations of Figure 3. Also, it may be noted that 
tha more pronounced decay [H] is .*emarkably similar to that of tha sodium 
decays of Figure 2. Clearly, there is some correlation between Na and H 
with the sodium decay obviously tracking that the H-atom. This behavior, 
similarly observed in a previous study of the laser induced chemistry of 
sodium, suggests the presence of a steady state distribution that readjusts 
as the radicals decay. 

It doe<> seem appsrent therefore that the initial rates of decay of sodium 
are not a measure of the progress of the three bodyNa+ 02 'tM reaction as 
suggested by Kaskan^ and McEwan axid Phillips^ and so cannot lead to a 
meaningful value for its rate ccmstant. Jf anything, the sodium decay is more 
closely related to the H-atom decay, which is predominantly controlled by 
the three-body H402*tM -• H02‘^M reaction in these flames. It seems probable 
tliat the observed dependence on O, and M concentrations cl these initial 
sodium rates of decay was a consequence of this H 4 O 2 +M reaction rather 
than an indication of the occurrence of the Na402‘<’M reaction. This further 
indicates the care required in interpreting such flame profiles and emphasizes 
the need to measure as many species in as wide a range of differing flame 
compositions and temperatures as is possible. 


[h], ( molecule cm"’) 
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Figure 4. H>atoxn concentrations calculated from the measured 
OH concentrations. 
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McEwan and FhiUipa'^ aasumad th» pr^cwunanca o£ N&O 2 and its 
aqoUibration via reaction (1) in tiMir brief study of Oiis system. Usii^ 
our more extensive data base we have tested this hypo^’l.esis by calculating 
the equilibrium constant of this reaction at points throi^hout each of the 
flames assuming that loss of sodium is due t<rtally toNa02 formation. The 
result is indicated in Figure 5. The expected straight line corresponding 
to a 0Q(Na»02) 55.9 kcal mol*^ is not confirmed. The departures in each 
flame clearly indicate that flame radicals are involved. A similar test of 
tiw dominance of NaOH via the balanced reaction (12) is illustrated in 
Figure 6 together with the expected equilibrium constant calculated from 
the JANAF tables. Behavior above and below 2000K appears to be significantly 
different. Above 2000K an approximate relationship does appear reasonable. 
However, below 2000K departures are such as to suggest the participation 
again of supplemental chemistry that involves the radical species whose 
ccmcentrations significantly vary with time. The system is obviously more 
complex than previously considered and.a more complete kinetic modeling 
is necessary. 

A model has been analyzed involving all the possible reactions of 

Na> NaO, Na02 and NaOH in an H 2 /O 2 flame environment. These are listed 

in Table n and schematically illustrated in Figure 7 . An examination of the 

possible reactions of NaH indicate that it will be negligible in oxygen -rich 

flames and is not included. Also, due to the very low concentrations of 

sodium, molecules containing more than one atom of sodium cannot form. 

No reliable value exis'':s for the ^iH£(Na02) and the effect of this variable on 

the respective reaction enthalpici? ‘s included. At first sight, the possibility 

of handling such a large number of reactions seems questionable due to the 

7 11 

fact that rate constants are available only for kj, k ^2 ^-12* * However, 

a closer examination immediately leads to various significant conclusions 
and indicates that more than half of the reactions listed play a negligible 
role. Comparison to the analogous reactions where H replaces Na shows 
that it is not unreasonable to expect reactions (3). (4), (7). (9). (10). (16) 
and (17) to have preexponential rate constant factors that correspond to 
hard sphere gas kinetic collision frequencies. Consequently, in the initial 

14 
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Table 71. Liat of p<^ntial reacUone of Na> NaO, 
Na02 aadNaOH in flames. 



Reaction 

AH®(2000K) kcskl mol” ^ 

1, .1 

Na+Oj'^J^ =Na02+M 

. (dJ + 1.16) 

2, -2 

Na02+H =Na+H02 

- (53.56 - dJ) 

3. -3 

Na02+H =NaO+OH 

- (45.98 - D°) 

4. -4 

Na02+0 =Na0+02 

- (61.54- D®) 

5, -5 

Na02+0H sNa0+H02 

+ (Dq - 9.78) 

6, -6 

Na02+H =Na0H+0 

- (64.51 - dJ) 

7, -7 

Na02+0H =Na0H+02 

- (80.07 - D°) 

8, -8 

NaOH+OH = NaO+H20 

+ 1.89 

9. -9 

NaOH+H =NaO+H 2 

+ 16,55 

10, -10 

NaOH+O =NaO+OH 

+18.53 

11, -11 

NaOH402 =Na0+H02 

+70.29 

12, -12 

Na+H20 = NaOH+H 

+41.90 

13, -13 

Na+OH+M = NaOH+M 

-81.23 

14, -14 

Na+H02 =NaOH+0 

-10.95 

15, -15 

Na+OH =NaO+H 

+43.78 

16, -16 

Na+H02 =NaO+OH 

+ 7.59 

17, -17 

Na-K52 =NaO+0 

+59.34 


Dq represents DQ(Na- 02 ) 


Reacticm nusoiber: + forward reaction, - reverse reaction. 
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Figure 7. Reaction scheme including all the possible reactions listed in Table II that couple 
Na, Na02> NaO andNaOH. Specific reactants involved are indicated on the 
connecting arrows. Reaction number refers to Table n. 
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aiMilysls all reaction ratet have been set at their gas kinetic equivalent 

val<M suitably modified in the ai^ropriate direction by the equiUbrium 

constant. The required equilibrium cmstant for Na02 has been calcula^d 

based on a 105 bent structure with r ^a-O) 0. 20 nm, rJO-O) 0. 13 nm and 

* -1 17 18 

vibratimial frequencies of 1094^ 391 and 333 cm ' Any lacertainties 
in this value are in fact (mly of specific concern to reactions (3), 

(4). (6) and (7) and should not introduce significant errors. 


An analysis of the time constants of these reactions confirms that 

the ms time franoe provides a sufficient time for the establishment of a 

steady state distribution. Closer examination of the kinetic model immediately 

indicates the implications of Husain and Plane's^ value of 1.0 X 10 cm^ 

“2 - 1 

molecule s for k|. Such a value predicts a severe depletion of Na that 
is too large to be compensated for unless DQ(Na- 02 ) is drastically reduced 
from its previously reported value of 55.9 kcal mol"^. Also the flux 
coupling NaOH to HaO is so large that their distribution will be described 
very closely by the balance of reaction (8), suggesting thatHaO will contribute 
at the most only a few percentage to the total molecular formation. 


A test of the participation of onlyNa/NaO^ andNa/HaOH reactions 
shows an unreasonable fit to the sodium curves of Figure Z, even when 
DQtHa-O^) is drastically reduced and the appropriate rate ccxxstants are 
varied within reasonable bounds. Clearly some direct coupling between 
Na02 and NaOH is required to adjust calculated to experimental N a values. 


Computer fitting of the model to the data has indicated that the most 
sensitive parameters are Dq ^a- 02 ) and k^. Unless the bond strength is 
reduced to 35 kcal mol , and appears fixed to better than + 5 kcal mol'^. 
and kj is reduced to 10 cm^ molecule ^ s it seems impossible to 
fit all flames to the data. Both shape and magnitude have to be reproduced 
which is quite demanding even though numerous parameters can be 
adjusted. 

Our preliminary optimum fit for three of the flames is indicated in 
Figure 8. Refinements are still in progress and we expect to complete 
this program very soon. 
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DOWSTREAM TIME, (ms) 


F4(ure 8. Optimized fit of the complete kinetic model to the experi- 
mental 8 odiurn concentrations. DS = 35 kcal mol~^, 
k| = 10-31 cm° molecule"^ 
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The calcttlatl<me indicate that of the varioua reacticma Na and NaO^ 
are coupled predominantly via reaction (1) and reaction (2) is negligible. 

Direct coupling between Na andNaO is negligiUle and reactions (15)— (17) 
can be ignored. Na and NaOH are conned by thm well established reaction (12) 
which is known to be efficient. Due to the fast coupling betwemNaO and NaOH 
through reactictt (8) the rates of (9)— (H) do n<^ have any c<mtrol in the model. 

Ih the other branches, reactions (3), (6) and (7) appear important and runs 
testily; the sensitivity to their rate constants are in progress. It seems un- 
likely that the model will be able to distinguish between the indirect con- 
version of NaO^ to NaOH viaNaO or the direct channels. 

It is clear that a simplified mechanism involving the seven reactions, 

(l)» (3), (4), (6), (7), (8) and (12), fits the data quite well. Large values for 
the rate constants of most of these a*re found necessary. 

Recently, the apparent discrepancy between Husain and Plane's 

Q 

measurement of and Kaskan's flame analysis has been realized by Jensen.^ 
Postulating thatNa02 plays a negligible role, he has suggested in a reinter- 
pretation of Kaskan's limited flame data that the reactions (12) and (13), namely 

Na+H20 NaOH+H 
Na+OH+M NaOH+M 

are sufficient to explain the main features of the flame data. A rather large 
rate constant of 5 x 10 T ^ cm^ molecule"^ s”^ is suggested for the three 
body formation of NaOH. This is two orders of magnitude faster than the 
corresponding H-analog reaction used in our analysis. Although open to possible 
uncertainty, it might be noted that also appears to be much faster than might 
have been expected. 

Testing Jensen's hypothesis on our larger set of data confirms our own 
conclusion that a more complete chemistry must be considered. Utilizing 
the rate constants suggested by Jensen shows a poor fit to the data, in both 
shape and magnitude, as indicated for two flames in Figure 9. Modifications 
to the rate constants of the two reactions, although improving the fit, shows 
an inability to provide an acceptable agreement. The reason for the inadequacy 
of the model may be illustrated from plots of the molecular distributions pre- 
dicted from our optimum fit to the more complete kinetic scheme. As seen 
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in Figure 10, although MaOH doe* coaalatently appear to be the dominaav 
molecular form aa euggeated by Jeaam, the role of Na 02 cannot be ignored 
particularly at temperaturea ci < 2000K. Conaequently, in ap^ of its weak 
molecular atructure, ita rapid rate of produetictt can maintain a significant 
ateady atate population. 

Further refinementa of our model, apecifically to examine the efifecta 
of significantly increasing the rate constant of as im^ed by Jensen 
and Jones^ now are in progress. 


original page is 
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